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Abstract 
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a,: 

Qh! The high energy behaviour of two-scale hard processes is investigated in the framework 

of small- £ models with running coupling, having the Airy diffusion model as prototype. We 
show that, in some intermediate high energy regime, the perturbative hard Pomeron exponent 



determines the energy dependence, and we prove that diffusion corrections have the form hinted 
at before in particular cases. We also discuss the breakdown of such regime at very large 



energies, and the onset of the non-perturbative Pomeron behaviour. 
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1 Introduction and Outline 

High energy hard scattering has received considerable attention in recent years. The essential 
problem is to determine the Green's function G(Y; k, ho) for gluon-gluon forward scattering, 
where k and k are the mass scales of the gluons and Y = \og(s/kk Q ) is the rapidity cor- 
responding to a center-of-mass energy squared s. The classic calculations done by Balitsky, 
Fadin, Kuraev and Lipatov (BFKL) U several years ago corresponds to an approximation (the 
leading series of powers in aY) where the QCD running coupling a s (t) is treated as a constant, 
a. In this case the high-energy behaviour of G is determined by the rightmost singularity in 
the ui plane, where ui is the variable conjugate to Y . This singularity at u s = ax™ = ^Xm is 
given in terms of the saddle point Xm of the function x{l) which gives the eigenvalues of the 
BFKL kernel. 

When higher order corrections [0, |3|] to the BFKL kernel are taken into account, the situation 
changes conceptually due to the running of the QCD coupling [§], |J. Running coupling effects 
mean that the saddle point of x is now a function of the scale t = \og(k 2 /A 2 ), u s (t) = a s (t)x m - 
Furthermore, u s (t) is not a point of singularity of the Green's function g w (t,to), although the 
value uj s (t) does control the F-dependence of G over a limited region of moderately large Y- 
values f5j. The rightmost singularity of g w (t,to) is at uj = up, is independent of t and to 
and determines the asymptotically large F-dependence of G, although for t and to sufficiently 
large this asymptotic behaviour will not set in until Y is quite large. The singularity at up is 
determined by non-perturbative physics. 

The fact that running coupling effects can change the character of the F-dependence of G 
is easy to see. In the fixed coupling limit the BFKL kernel gives a contribution proportional 
to aY each time it acts. In the running coupling case, the contribution is proportional to 
a s (t')Y ~ a s (t)Y + b\og(t' /t)a 2 (t)Y , when expressed in terms of the fixed external scale t of 
the scattering. However, the contribution of a single running coupling term vanishes, because 
the average value of log(t'/t) is zero since the probabilities for t' > t and for t < t are equal in 
fixed coupling BFKL evolution. At the level of two running coupling contributions, one gets 
the average of a^(t)Y 2 log 2 (t ft') ~ a^Y 3 |B|, [7], |^, M, and this contribution exponentiates. This 
simple perturbative argument is valid so long as a^Y 3 <^C 1, but is difficult to extend to values 
of Y > £ 5 / 3 . And it is exactly in the region £ 5 / 3 < Y < t 2 where the most dramatic running 
coupling effects on BFKL evolution take place. 

In the present paper we calculate, starting from some small-x models, the diffusion and 
running coupling corrections to the hard-Pomeron exponent. Basically, we prove the validity of 
the leading running coupling corrections hinted at before 0-0, in the full range t <^Y <^t 2 , 
and we discuss some features of the regime Y ?Z t 2 . 

Because of the conceptual complexity caused by the running of the coupling, it is useful to 
have a simple model where the essence of running coupling effects are present and yet a rather 
explicit discussion of the F-dependence, in the various regimes of Y, can be given. Such a 
model was introduced some time ago |J, which takes into account the running of the coupling 
as well as diffusion in momentum scales in terms of the quadratic dependence of x{l) about 



the saddle point at 7 = 1/2 (Sec. 2). Here we study the perturbative behaviour of this (Airy) 
model and of a simple generalization where x{l) is represented as a sum of two poles in 7 0. 
These two models give identical results for G in the region t < F < t 2 , and we believe they 
represent QCD accurately in this region, as outlined below. 

The solution of the Airy model is given in terms of Airy functions for the perturbative part 
of the evolution and in terms of a reflection coefficient, determined by the way the running 
coupling is regularized in the infrared, for the non-perturbative part (Sec. 2). The singularity 
at up resides in the reflection coefficient S(u) 0. 

When t <^Y <^i t 5//3 the behaviour of the perturbative part, Gp(Y; t, to), of G, is determined 
by a saddle point of the cu-integral at u s (t) giving Gp ~ exp(u s (t)Y). In this region of Y running 
coupling effects play a minor role. When Y reaches t 5 ' 3 the saddle point at u s (t) no longer gives 
the dominant contribution to Gp. Nevertheless, the largest term in the exponent governing the 
F-dependence of Gp is still u s ((t + t )/2)Y, with the next largest corrections being given by the 
"diffusion" term and by the Y 3 /t 5 term discussed above (Sees. 3, 4). This result is summarized 



by Eq. (^1|) and by Eq. ([50]) of the text, which give the dominant behaviour throughout the 
region i < 7 < t 2 . When Y reaches t 2 the character of the solution changes drastically. For 
Y ^> t 2 the magnitude of Gp goes as exp(u\/Y) where uo is a pure, t- independent, number (Sec. 
5). This behaviour comes from two complex-valued saddle points having uj ~ 1/y/Y, whose 
exact magnitude is model dependent (Sec. 6). Because the saddle points are complex there is 
an oscillating prefactor in Gp, which cannot be given a sensible physical interpretation, and 
calls for non-perturbative contributions to take over. 

The physical interpretation of our main results, Eqs. (21]) and (]50|), seems clear. For 
i<F< £ 5 / 3 , Gp behaves - with good approximation - as in the fixed coupling case, having a 
magnitude proportional to exp(u s (t)Y), with a spread in t — to = At given by (At) 2 ~ u> s (t)Y. 
At Y ~ t 5//3 running coupling effects become more important. For t fixed, Gp(Y;t, to) reaches 
a maximum value proportional to exp(u; s (t) + i] 3 /3) with 77 ~ Y/t 5 ^ 3 , while At is well fixed at 



a value ~ Y 2 /t 3 , with only small fluctuations of size yY/t from that value. Thus the values of 
At are no longer given by pure diffusion, but now (for a fixed t) to is being pulled in the infrared 
where the coupling is large. When Y gets as large as t 2 the preferred value of At becomes as 
large as t and the perturbative part of the model ceases to make physical sense. 

Finally, a comment on the accuracy to which we have calculated Gp. It is convenient to 
write logGp = uo s Y f(Y 2 /t 4 ,At/t) + small terms. We have calculated terms in / of size 1, 
Y 2 /t 4 , At/t, and (At) 2 t 2 /Y 2 , as given in Eq. [H] In the dominant region of At ~ Y 2 /t 3 all 
these corrections are of size Y 2 /t 4 . We certainly expect further corrections to / of size Y 4 /t 8 , 
as discussed in Sees. 3 and 4. So long as Y <C t 2 we believe we have the dominant terms 
in the exponent and so have the essence of the growth of Gp with Y and its dependence on 
At. However, we only have some preliminary ideas on how to calculate corrections beyond this 
region (Sec. 5), and thus we likely do not have all the large terms in the exponent of Gp. 



2 The gluon density in small- £ models 

We consider in this paper a particular form of the 2-scale gluon Green's function which has 
been established for the Airy model and for the truncated BFKL models ||. By defining 

G(Y;k,k ) = -±- J ^ Y g u (t,t ) , (1) 

k 2 

where Y = log ^J- and tj = log ■$, y w takes a factorized form for t > t , namely 

g u (t, to) = F u (t) fe(to) + 5(w)^(t )) , (t > t ) , (2) 

where the various terms are defined as follows: 

- -^L(t) is the "regular" solution of the homogeneous small-x equation being considered, which 
vanishes exponentially for large t values; 

- ^L(to) is an "irregular" solution, which is instead exponentially increasing with to; 

- S(u>) is a "reflection" coefficient, which has been explicitly constructed in some cases ||, |J 
and is dependent on the strong coupling region, e.g., on how the t = Landau pole is smoothed 
out or cutoff. 

While the explicit form and size of the non-perturbative S(u) part is dependent on the model 

- in particular on the number of poles taken into account in the effective eigenvalue function 
[Of, the perturbative term is unambiguously defined in the large-t region, and is supposed to 



be calculable in a realistic small-x framework [ lQ| . For this reason, most of our analysis will 
concern the perturbative term. 

In order to understand better the meaning of Eq. (|2|), let us derive it explicitly in the Airy 
model. The defining equation for the Green's function is, in operator notation, 

(u - a,K)g u = 1 (3) 

where K is in general an integral kernel, and 6t s (t) = Nccx s (t)/n = 1/bt is the running coupling. 
The Airy model obtains by assuming that K is scale invariant, and described by a quadratic 
eigenvalue function 

x (^)~ Xm + - x ^( 7 --) 2 , ( Xm ,x'L>o), (4) 

which is a reasonable approximation around the minimum of Xi which is taken to be at 7 = |. 
Here 7 is a variable conjugated to t by Fourier transform. Therefore, by using Eq. ([|) in 
t-space, Eq. (0) becomes 



uo - uo s (t) (1 + Dd 2 t 



g^t,t ) = 6(t-t ) , (5) 



where u s (t) = Xm&s(t) is the hard Pomeron exponent, and D = ^Xm/Xm ls the diffusion 
coefficient. 



In other words, the Airy Green's function satisfies a second order differential equation in 
the t variable and has, therefore, the well-known form 

9u,(t,t )=to(r*(t)rt;(t ) e(t - 1 ) + (t <-► t )) (6) 

where 3 rR (t)(toJ-' l , j(t )) is the regular solution of the homogeneous equation in (||) for t — ► oo 
(of the adjoint homogeneous equation for t — > — oo). Equation (|j) is the basis for Eq. (0), but 
should be better specified by smoothing out or cutting off a s (t) in a region t ^ i , where t > 
defines the boundary of the perturbative regime a s (t) = 1/bt. Depending on such procedure, 
the left-regular solution can be evaluated for large to values in the form 

*S(*o) = ^(to) + S(u)^(t ) , (7) 

where T 1 is irregular for to — > oo, and S{u>) is a well-defined reflection coefficient. We have 
thus derived Eq. (0) with T = T R and T = t§T l . A similar derivation holds for the truncated 
BFKL models with n poles, and in particular for the 2-pole model (Sec. 6). 

In the Airy model, Eqs. (|6|) and (0) can be given in explicit form. In the perturbative region 
t > t we can set, by Eq. (J5|), 

and furthermore, 

W V X m / Az(£) 

where Az(£) is the Airy function, 

5i(0 = e~,4z(£e— ) + e~ Ai(^e^) (10) 

is the irregular Airy solution, 



1 / ,\ 2 



'-©■(-^-^(^rfesH (n) 

is the Airy variable, and £ is its value for t = t, at the boundary of the perturbative region. 
Note that, given the delta-function source in (H), the regular and irregular solutions in Eqs. flD 
and (H) must have a well defined Wronskian. 

From the explicit expressions fl8D-(jl0"D it follows that J- w and J- u are analytic functions of u, 
showing an essential singularity at u = only. Instead, the reflection coefficient S(u) - quoted 
in Eq. (g) in the case a s (t) is cutoff at t = t - shows uj singularities in Recu > at the zeros 
of the Airy function, meaning that the leading Pomeron singularity Up actually occurs in the 
non-perturbative term in Eq. (§). However, the latter is suppressesd, at large t , by the ratio 
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Figure 1: Integration contour Cr + C s for the Green's function decomposition in Eq. (P"2"|). 
Cr{C s ) refers to the Pomeron (perturbative) contribution. 

Ai(C,o) / Bi(C,o) , meaning that the perturbative term may actually be more important at large 
scales and intermediate energies. 

By using the decomposition in Eq. (|2]) we can rewrite Eq. ([!]) as a sum of two terms 

G(Y;k,k ) = G P + G R , (12) 

where 

G R (Y; k, k ) = ^- f ^-.e" Y S{u) F u (t)F u t ) (13) 

rhrhQ JCji ZiTX% 

carries the (power behaved) Regge contributions (Fig. 1), the leading one being the (non- 
perturbative) Pomeron, while 



G P (Y;k,k ) 



1 



du 



kko JCs 2iri 
1 

kko 



e uy Tu,(t)^(t ) 



exp(uj s (kk )Y + 



(14) 



corresponds to the "background integral" and is characterized by the two-scale exponent u s 
that will turn out to be determined by a s (kk ) (Section 3). Our goal in the following is to 
analyze in more detail the ^-dependence in Eq. (0), by determining the regime in which the 
exponent uo s is relevant, and the magnitude and form of diffusion corrections to it. 



3 Diffusion features of the Airy model: a heuristic 
approach 



We have just clarified that the two terms in the decomposition (Q) generate the F-dependence 
in G(Y; k, kn) in a quite different way. The non-perturbative (Pomeron) term provides just a 
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Figure 2: Contours in the v = 7 — 1/2 plane for (Cr) the regular and (CV) the irregular solution 
in the case of (a) the Airy model and (b) the two-pole collinear model. 



Regge-pole behaviour ~ exp(upY), while the perturbative one is analogous to a "background 
integral" (Fig. 1) and will generate a non-trivial exponent only if the small-cu oscillations of T 
and J- are kept in phase by the cu-integration. By writing, for the Airy model, 



g P (Y;t,t ) = kk G P 



du wy vrt {2bu\3 

^— e — ( — ) MO B <Za) , 



v ^Xl 



(15) 



we expect, for At = t — to "C t that phase relations are kept only if £ and £0 are kept finite for 
large Y. By the definition (|TTD, this implies that 10 — uj s {t) ~ uj — u s (t Q ) = 0(t~i) are small 
parameters. Furthermore, in this region, by Eq. (|Tl~D, 



£~ irM 



uj 



Us{t) 



(16) 



is just linear in u — uj s . By replacing the linearized expression (|16l ) into Eq. (|15|) , we can rewrite 
it in the form 



g P (Y; t, to) ~ n(Dt)-* e^ Y / -1 e^4i(£)£i(£ - 5) 

where we have introduced the parameters 

77 = £>3 y u s (t) rt ~ yr§ , 5 = At{Dt)~§ . 



;i7) 



The integral in fll7|) can now be evaluated in a heuristic way by introducing the integral 
representations (Fig. 2a) 
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(19) 



lc++c- 

and by performing the ^-integration under the integral, to yield a delta-function. We thus 
obtain 

n 3- 



1(77, 5) = f P, e*» Ai(0 Bi{Z - 5) = [ + ™ p± exp Lv\ + ( V 2 + 5)v A + V S + 



V 



By inserting Eqs 

9p(X,t,t 



2ityfWrj 
0|) and flTJ|) into Eq. (|I7|), we finally get 
1 



IOC Z71 % 

1 (if r]S 5 2 \ 
exp ( — + — — ' 



exp 



2 4r?y 

(At) 2 



W ' ( * )y < 1 + ¥)-40*(i)Y' 12 



r? 3 



(20) 



(21) 



/47rDw a (t)Y 
which provides the diffusion and running coupling corrections to the hard Pomeron behaviour. 

There are two features worth noting in Eq. (^l|), in comparison with the customary expres- 
sion with frozen coupling. Firstly, the exponent uo s is corrected by a term linear in At = t — to 
which provides the symmetrical argument uo s (J^-) in the running coupling, as is appropriate 
for the factorized scale ~ kko already introduced in eq. flT2|). Furthermore, the exponent carries 
the diffusion correction rf ~ Y 3 /t 5 , which is of relative order Y 2 /t 4 compared to the leading 
term u s (t)Y. This correction was obtained as running coupling effect in ||, [pi, Hand confirmed 
in the models considered here under the assumption 77 <C 1, or Y <C ts. 

The question then arises of the boundary of validity of the heuristic argument presented 
above. We shall see in Section 4 that the assumption 77 <C 1 can be relaxed, and replaced by 
77 <C £3, or Y <C t 2 . For Y ^ t 2 on the other hand, the linearization in Eq. ( |I4"D breaks down, 
and the integral in Eq. (13) enters a new regime in which it first decreases, and then starts 
oscillating (Section 5), so that the phase relations are lost. 



A first hint at such behaviour is obtained by replacing the ansatz in Eq. ( |2T| ) in the diffusion 
equation 

g P (Y,t,t )~uj s (t) (l + Dd?) g P (X,t,to)- (22) 



BY 
In fact, the ^-derivative of the exponent E(Y, t, At) = log gp is given by 



dE(Y,t,At) 



1 



-^+S(t)[l + -) 



At\ 



(At) 2 _8_ 
WuoY 2 + dY 



dY 2Y 

and should match the right-hand side of Eq. ([22]), which is given by 



rf 

Y2 



(23) 



u B (t) 

dE 
~dt~~ 



1 + D 



' 9 /H 2 + Q2E 



At 



2Dlo s Y 

7 



t ' dt K w 



1 

2£ 



(24) 



If we keep terms of relative order t/Y, At/t, Y 2 /t 4 (such terms are all of the same order for 
fixed values oirj), we find that Eqs. (53) and fl24|) are indeed consistent, provided 



dY Y1~ D 4*2 ' [Zb) 



thus reproducing the expression (|18D for 7](Y,t). But, because of (|25|) , we also generate from 
Eq. (|22| ) terms of relative order Y 4 /t 8 which are subleading only if Y <C t 2 , and do not check 
with Eq. (|23|) . Therefore, the heuristic argument breaks down for F = 0(t 2 ). 



4 Detailed analysis of the regime t <Y <t^t 2 

In this section we give a more complete evaluation of the integral (|i~5|) which defines gp and in 
doing so we confirm the result (|21|) and its validity boundary, obtained in a heuristic way in 
the previous section. 



4.1 Choice of integration contour 

We begin with the cu-integration contour in (|i~3"D being parallel to the imaginary u;-axis with 
Reu = u) s (Fig. 3). In order to effectively separate leading from non- leading behaviours in ( |T3"D 
it is convenient to use different forms of the product Ai{£) -Bi(£o) for Imtu > and for Ivauo < 0. 
To that end we write 



M0B i ^)=T 2 -^K^I[ 



^(pe-)-i^(p) 



(26) 



where the upper (lower) sign will be the form used in (^) when Imu; > 0(< 0). Equation (|26] 
follows from (|8|) and ( |T0D along with 



Thus 

with 



K u (e™() + K u {e~™Q = 2 cos nuK u (C) ■ (27) 

g P (Yt, t ) = g P (Y; t.t ) + R(Y; t, t ) (28) 



UriY-.t-t,) = J T J^- /( " _i)( " ~^ <f*K k {\&) KiMe^) (29) 



and 



*™« = I ± ^d ^-v^- 1] e " YK ^ ^ < 30 > 



where the upper (lower) sign is to be used when Imu; > 0(< 0). 
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We shall first show that R can be chosen small compared to go by a judicious contour 
deformation. Since the integrand of ( f2"4| ) decreases for positive real £, we are led to deform C s 
around the real £ axis to reach a point £1 > 1 to be defined below. In estimating the size of R 
we shall take £o ~ £ so that the full exponent appearing in ([!(]), in the region where |£| 3> 1 
where the asymptotic form of K\ can be used, is 

8 = u s (t)Y + {u - u s (t)Y - ~Z% . (31) 

2 

In the region where £ <C ts 

£-u s Y~tr)- ^i . (32) 



Thus we can make the second term on the right-hand side of (p2|) dominate the first term at 

£ = 6 if 

6/V » 1 • (33) 

We note that (|33]) is possible, keeping £i <C t 2//3 , so long as rj <C frs or, equivalently, so long as 



7 < i 2 . We anticipate gp being of size e^ sY so that if we are able to choose an integration 
path in ( f24"|) having Re|^2 ^> (u — u s )Y, then R can be neglected. A contour of this kind is 
C*2 in Fig. 3. It is basically a deformation of the contour C s at Keu = uo s in order to have 
£i as starting point and to depart from it with Arg(£ — £i) < §■ Its basic property is to lie 
completely within the regions |Arg£| < |, such that the product of regular Airy functions is 
damped (Re£a > 0). It is then easy to convince oneself that Re(|£a) 3> Re(u — io s )Y on a 
contour of this kind, so that the contribution R in Eq ( |3TiD can be neglected compared to gp so 
long as Y >C t 2 . 



4.2 From the cj-integral to the ^-integral 

Our task is now to evaluate ( pS| ) and to specify once again a convenient path of integration in 
£, that will turn out to be different from that chosen before. We are able to evaluate ( p9|) only 
in a linear approximation of £ and u — u s (t), and we first turn to a more complete justification 
of this approximation. When |£ | and |£| are large, the terms appearing in the exponent in (|29|) 
are 

£ = uY-^ + p . (34) 



This can be written as 



uj \ 3 2 3 2 , t _ / u; \ 3 , 3 



e= w . W y + e,l-j -jf + jtt-^ -))>, (35) 



with 5 as given in (|18D. The linear approximation gives 



£ lin = w,(t)y + ^ - ~ £§ + !<£ - 8)1 , (36) 



9 




Figure 3: Integration contours in the £ plane, cut at £ c ~ t 2 ^ 3 . C s corresponds to Keco = uj s . 
The dashed (dotted) contour C\ [C-i) corresponds to the deformed contour defined in the text. 
The corresponding ones in Im£ < are not shown 
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In order to see how close £ and £\[ n are, we deform C s to the contour C\, in this case (Fig. 
3). The latter, departing from £1 and reaching the origin, is chosen in such a way as to run 
over the imaginary £ axis in the whole linearization region |Im£| < ts, and to approach C s 
later on in the Reu; < u; s region. It is straightforward to see that the regions of C\ on the 
real and imaginary axis are the dominant ones in the integral. In fact, from Eq. ( J36l) one sees 
that Re(£ lin — uj s Y) vanishes at £ = and decreases steadily below zero for £ on the imaginary 
axis. Furthermore, in the asymptotic region where C\ starts approaching C s the Airy phase is 
negligible and the integrand in Eq. ( ^9|) is damped because of Keu < io s . 

Now we want to replace £ by £n n in the exponent. We can estimate the error in the exponent 
by comparing (|35| ) and fl36|) in the region < £ < £i. The maximum deviation occurs for £ = £i 
and the size of the deviation is 



AS = £ - £ 



lim 



£iv 



& 



- 1 



(37) 



From (O) one sees 



- 1 



and 



- 1 



are of size £it 2//3 , so that A£ has terms of 



size £ 2 ??t _2 / 3 and £■/ At ft. Taking £i = Nrf for N large, the smallest we are allowed to take 

£i, and taking At ~ (Dt)3?f as region of At where gp, as given in (0), take its maximum 

value we find 

5 y 4 

A£^X, , "UsY-jr . (38) 



rf 

2 



Thus we can expect our final result and have errors of size A£ in the exponent. 



In changing the exponent £ in (|29f ) to £u n we take the same contour as indicated in Fig. 
3, but now the integration is continued along the imaginary £-axis up to £ = zoo. Using the 
linearized form for d£ also induces an error in the prefactor of the dominant exponential of size 

2 



OC 



t 4 



compared to 1. Thus, finally, the integration which we need to evaluate is 



g P (Y;t,t ) = T [Z)t]-3e ±i ie' Js(t)y / ci£ e^ Ai(g)Ai((£ - 5)e T2i % 



(39) 



where we have used (g) to express the K\ functions appearing in ( p3|) by Airy functions for 
reasons which will become clear in a moment. The integration is taken with the contour C\ 
extended to ±ioo, with, as usual, the upper (lower) sign referring to Im£ greater (less) than 
zero. 



4.3 Evaluation of the £-integrals 



In order to evaluate the integral in ([39]) it is convenient to define 

4 = / di e*P ± (£) , J± = / d£ e*Q±(£) 



(40) 
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p±(0 = u(t)v ± {t - 6) , g ± (o = £zU{t)v ± {t - 5) 



with 



(41) 



aT 2if 



U(0 = Ai(0 , \4(0 = Te ±l 3 Ai((£ - 5) e 
Using ^yC/(0 = */"(£) = £U(£) and tf£(f - 8) = (£ - <*)Kt(£ - 8) one can show that 

p'" = [4^P - SP]' -2P- 2SQ , Q" = 2(£P)' - P - SQ 
from which, by integration by parts, one finds 

21 = -r 1 [A^--5-r 1 2 ]I-2SJ+e r >^P2^i,v) , 



J{6 + r) 2 )+I 



81 
-2/7— + e^Q 2 (M 



where the ± indices have been suppressed in (^)-(45). Also 

Pi = p"{Zi)-vp'{Zi)+v 2 p{Zi)-(4Zi-6)p{Zi) 

Q 2 = Q'^ 1 )-2^P(^)- V Q^ 1 ). 
Using Eq. (H), one easily finds 



d 1 
V *j + 2 - 



(5 + r] 



2N2 



The solution to (HB) is 



'± 



>fn 



exp 



4r/ 



2n2 



'± 



At] 2 



(S + r ] 2 )P 2± -26Q 2± 



S. 



± 



?7 3 (5 - rf) 



4r/ 



C± + 



^ drj 



exp 



r/' 3 (5-i]' 2 ) 



2\2 



4?7' 



"70 ,/l] 

with C-t a constant to be determined. Referring back to ([39]) one can write 

g P (Y;t,t ) = {D 2 tyU^)y (/+-/_). 



s±fa' 



(42) 
(43) 



(44) 



(45) 



(46) 



(47) 



(4* 



Now the combination S+ — S- occurring in the integrand of the particular solution of (|4"7 ) 
which enters in fl4*8|) always involves the combination 



^i(~(6 " *)* e-) + K^fo - 5)1 e + n = X*(~(6 - 5)1) oc exp[-|(& - *)*] 

so that the particular solution can be ignored. Although we have focused here on rj large the 
procedure also should be valid when 97 <^C 1. For rj -C 1, the first factor in (|47|) , -^= exp[— |-], 
gives the usual BFKL diffusion which means that C + — C_ can be determined by evaluating 
I + — I_ when 5 = and when rj is small. One finds 



C+-C- 



(49) 



12 



giving exactly the result (21) 



///.(>•://„) = , = e ".W-fts-W+4 (50) 

AnDuj s {t)Y 

In arriving at this result we have chosen to expand the integrand of (p9|) about £ and t 
and then integrate over £. Alternatively, if the expansion had been about £o and t with the 
integration over £ the result 

g P (Y; t, t ) = , 1 e ».to)r-4sl*>-<ir+4 (51) 

kirDu a {to)Y 

would have been obtained with r/ = Yuj s (t )D3t 3 and with <5 = (t — £) (Dt ) 3. The 
difference of these exponents is 

5 r) -8r] 1 (5 2 S 2 \ | ?? 3 - ry 3 _ A 



2 4 V r/ W 12 

When f < 1, 

Using 5 oc ?? 2 as the important values of 5 in (0), one easily finds 

5 y4 

A«> ^(-rr) ( 54 ) 

J3 t 

which is exactly the size of the error estimated in fl32|). This conclusion is consistent with the 
estimate -^ = 0(^s") obtained from the heuristic argument based on Eq. (p2j). 



5 The very large Y regime 

In the previous section we have exploited the decomposition (|2lf ) of the Green's function in two 
terms, one (Eq. fl29|)) with product of Airy functions out of phase (~ K '(£) K \e T ™ £ )) , and the 
other in phase (~ K(£)1T(£ ), Eq. (|30|) ), corresponding to the exponents 

£(u,,Y)=u;Y-~^±~£f . (55) 



We have shown that, for Y <C t 2 , the integral is dominated by Reu ~ u s , the contribution ( |30| ) 
being negligible on a contour with a slight deformation in Recu > cu s . On the other hand, for 
Y ^> t 2 , it becomes profitable to distort the contour in Recu < u s , where the contribution 

3 

becomes actually dominant, because Re£a may take negative values. 
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In fact, for small enough u in Iinu^O, the exponent with (— ) sign in Eq. (|55|) takes the 
form 



£+~ujY ± — 



U tu s (t) 



3 Duo 



with derivatives 



d£± U Xr, 



d 2 £, 



± 0» Xr. 



du ' ' 3 fev^cu 2 ' du> 2 3 fev 7 ^ 2 

Therefore, there are saddle points of the exponent (Ef|) at 



(56) 
(57) 

(58) 

, / ( 59 ) 

Since (d 2 £±/duj±)~2 <C w± the gaussian integration, in (|24"D, about the saddle points should be 
an accurate evaluation of the integral. One finds 



with 



u± 



d 2 £ ± 
du) 2 



3 V ; UyVd 



, ., '3by/D\h * 

1=F«) ( — ) Y* . 



gp(Y;t,t c 



to f 36 
3v^ UxmDW,/ 



3 

\ 4 



cos(o;vy 



I) 



jVy 



(60) 



with u) 



-,2 



8Xm 



It then appears that the perturbative solution becomes smaller than predicted by the ex- 
ponent uj s (kk )Y (because VY <C uj s Y for Y ^> t 2 ) and furthermore it starts oscillating, thus 
loosing positivity. It is no surprise that, at such large Y values, our perturbative calculation no 
longer has physical sense. After all, when Y approaches t 2 from below it is clear from ([H]), or 
(|50|), that At approaches t at which point one expects the singular potential evident in (|5|) to 
become troublesome. That is, when Y approaches t 2 the cutoff at t is clearly necessary. Two 
questions then arise: firstly, whether we can still describe the behaviour of gp in the intermedi- 
ate region F oc t 2 ; and, secondly, at which Y values does the non-perturbative Pomeron part 
really take over. 

We are unable to answer either question in detail. However, for the pure Airy model 
described by Eq. (|22|), we can provide a partial resummation of the corrections to the exponent 
of relative order (Y/t 2 ) 2n as follows. Referring to Eq. (24) we first neglect, in the large Y 
regime, the term -^ compared to (^) 2 ; then we consider the exponent E(Y; t, At) around its 
maximum in At so as to neglect its At derivative, and we take the ansatz 



E(Y; t, At n 



u a {t)Y f(z) 



D X Y - = ^- 

t 4 w,r 



(61) 



which is supposed to describe the Y 2 /t 4 dependence. Finally, by replacing Eq. fl6"T|) into Eq. 



(22), we get the non-linear differential equation 

f(z) + 2zf{z) = l + z[f(z)+4zf(z)} 2 



(62) 
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which is a sort of Hamilton- Jacobi limit of the diffusion equation 



For z < 1 (Y < £ 2 ), the iterative solution to Eq. (H) is f(z) = 1 + z/3 + 2z 2 /3 + ..., which 
yields the if /3 term of Eq. ( |50| ) and the first non-trivial correction to it. On the other hand, for 
z ^> 1 (Y ^> t 2 ), Eq. ( |6^) still makes sense, with a solution f(z) ~ 2~ 1//4 + O^ -1 / 2 ), in order 
to compensate the term 1 in the r.h.s. Therefore, since z ex Y 2 /t A , the exponent in Eq ( joT] ) 
turns out to be in agreement with the saddle point estimate in Eq. fl60|). Note, however that 
while the behaviour Q5DD is supposed to be universal for Y <C £ 2 (Cf. Sec. 6), the resummation 
in Eq. fl6"2|) and the large Y regime in Eq. ( |60|) are typical of the Airy model, because they 
involve large £ properties of Eq. (^) . 

We wish now to compare the large Y behaviour of the perturbative term in Eqs. (]2l|) and 
fl50|) - with its intricate transition to the behaviour ( p0|) just discussed - to the non-perturbative 
Pomeron term. The latter can only be discussed in a model dependent way. In the Airy model 
defined in Eqs (§) and (|9|) it is straightforward to show that 

gR iY- tM ) = e ^^£(^) 2/3 A - (fr) a exp[wp y _ W^V' 2 ,, (63 ) 

WP V Xm > 3 \XmDJ 

where £p ~ tipup/xmD) 1 ^ and we have assumed c<jp ^> cu s (t), for definiteness. It is clear from 
Eq. ( |63|) that the Pomeron term is exponentially suppressed with respect to the perturbative 
one for large t. However, it takes over at very large energies, such that 

4t 3 / 2 / hup \ V2 
that is, even before the region Y ~ t 2 is reached, as already pointed out for collinear models in 

The above estimate of Y c is model dependent in several ways. In fact the weight of the 
Pomeron term depends on the small- a; equation being used (cf. Sec. (6)) and, within the given 
equation, on the way the strong coupling region is smoothed out or cutoff [|Kj. Furthermore, 



unitarity corrections may affect it, and there is no consensus on how to incorporate them. 
However, the basic qualitative feature underlying all models is that "tunneling" to the Pomeron 
behaviour occurs at some Y c < t 2 , i. e. before the perturbative calculation looses sense, thus 
insuring cross-section positivity. 



6 Extension to collinear small-x models 

The above evaluation of perturbative diffusion properties of the Airy model (Eqs. ( FJip and (150)) 
can be generalized to other small-x models, for which the expression ( |T2D of the perturbative 
Green's function is sufficiently explicit. 

For instance, a simple two-pole collinear model || is provided by the effective eigenvalue 

function 

11 

X(7, uj) = - + — , (65) 

7 1 + u — 7 
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where we have introduced the upshift [12], but no further subleading terms. The corresponding 
Green's function is defined by [DJ 



[w - X(j + 5 *' w )^] &>(*> *o) = <J(t - t ) 



(66) 



which basically provides a second order differential equation for g u . Therefore, the solution is 
of type in Eq. @, which, for t > t > i, reads 



9uj(t,t 



- S(t - to) + ^t— Wj, i 



'l + u/)t 



x 



W_±_ n 

bu'2 



-(l + w)t ) + 5(cj)iy^ i ((i + u%) 



(67) 



Here Wk tm (z) is the Whittaker function, satisfying the Coulomb-like differential equation 



d 2 k m? 



dz 2 



+ 



+ - 

z 2 4 



W hm (z) = 



(68) 



and its "irregular" counterpart "W_fc jm (z)" is more precisely defined by the real analytic con- 
tinuation 

W- k {-z) = - (e ink W. k {e m z) + e- ink W_ k {e- i77 z)) . (69) 



Such solutions admit also the 7-reresentation JlO] 



±(l+w)t 



W:J_ I/O, — 



d'j 



e (7-|*-5b X (^) 



(70) 



where Cr(Cj) is the contour for the regular (irregular) solution shown in Fig. 2b, and xili w ) 
a_ 

97" 



d -X(l,cu) 



The energy dependence of the perturbative function Gp(Y;t,t ) in Eq. ( |i"2|) was already 
investigated in Refs. JE(|,[ll| and is characterized by various regimes. 



1. In the collinear limit t ^> to > 1, and Y ^> 1, Gp(Y; t, to) is dominated by the customary 
anomalous dimension saddle point, which exists in the energy region 



where u s (t) 



(l+uj s )bt ~ bt 



^|^ ~ u > uM » u.{t) , (Y = log ±) 
is the saddle point exponent. 



(71) 



2. In the diffusion region At = t — to <S t, and 7 > i > 1, the important cu values drift 
towards u s (t) ~ co> s (t ), and the asymptotic behaviour of the Green's function matches a 
properly chosen Airy model. 
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In fact, if t and - are both large, but uj — u s (t) — uj — uj s (t ) is small, the phase in Eq. (|70| ) 
is finite only if (7 — ^) ~ (buj)s is small also. This means that we are probing the region close 
to the minimum of x{l)-> m such a way that 



1 . / ,\ . 2 



is finite. This is precisely the region where the Ws become asymptotically Airy functions, with 

parameters 

4 y" 4 4 

1 + ^s 2y m (1 + cj s ) 2 W 

For this reason, the linearized expression flT5| ) holds in the present case also, under the same 
conditions. 

The above argument shows that, in the region t <3^Y <^t 2 , and (t — to) €t, the diffusion 
corrections of Eq. ([19|) are valid for the two-pole model also, and for any truncated model 
where the 7-representation ( |70[ ) holds for both the regular and the irregular solution. 

Of course, for Y ^> t 2 , the important uj values become much smaller than uj s (t), and the 
phases in Eq. (|70|) became stationary at y(7,c<j) = 0. For the Airy model this occurs at 
v = 7 — I ~ ^T^' anc ^ P rov ides the asymptotic behaviour used in Eq. flSfip. For the two-pole 
model, instead, 7 drifts to ±zoo, and this provides the asymptotic formulas 

(2m 
g^ ~ exp I ± — 



g P (Y;t,t ) ~ expwv 7 ! 7 cos^x/F + -) , (74) 

8 

with a) = -y. Although the precise exponent is different, the qualitative behaviour ( [74] ) is the 
same as for the Airy model. This follows from the analogous structure of the 7-representation 
with its exp( £2S ^) behaviour. 

In conclusion, in the range i ^ F < t 2 , the perturbative calculation shows a universal 



behaviour, described by Eq. fl2T|) or Eq. (50 ) , where the only model dependence lies in the 
parameters y m and D, describing the hard Pomeron and the diffusion coefficient. For Y ^ t 2 
on the other hand, the perturbative behaviour is more model dependent and finally starts 
oscillating, thus loosing physical sense. 

At some intermediate value 7 = 7 C (t < 7 C < t 2 ), the non-perturbative Pomeron takes 
over. We do not have a reliable model for such transition. Is it an abrupt "tunneling" effect 
0, as in the two-pole model (with up 3> ui s (t)), or is it instead mediated by a long diffusion 
regime, as perhaps expected in unitarized models {up <^ 1)1 This is an open question, which 
deserves further investigation. 
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